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Abstract. - Cavities with periodically oscillating mirrors have been predicted to excite photon 
pairs out of the quantum vacuum in a process known as the Dynamical Casimir effect. Here we 
propose and analyse an experimental layout that can provide an efficient modulation of the effective 
optical length of a cavity mode in the near-infrared spectral region. An analytical model of the 
dynamical Casimir emission is developed and compared to the predictions of a direct numerical 
solution of Maxwell's equations in real time. A sizeable intensity of dynamical Casimir emission is 
anticipated for realistic operating parameters. In the presence of an external coherent seed beam, 
we predict amplification of the seed beam and the appearance of an additional phase-conjugate 
beam as a consequence of stimulated dynamical Casimir processes. 
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' Introduction. — On the basis of quantum field the- 
ory, a neutral planar mirror moving with a non-uniform 
deceleration is expected to convert the zero point fluctua- 
tions of the quantum vacuum of the electromagnetic field 
itito real propagating photons [T]. This so-called Dynam- 
ical Casimir Effect (DCE), predicted more than 30 years 
clgo, has proved extremely hard to experimentally observe 
due to the extreme weakness of the emitted photon flux 
fbr realistic configurations. Resonant enhancement of the 
DCE intensity inside an optical cavity with oscillating mir- 
rors was investigated HHi], still a sizable DCE emission 
requires a relativistic motion of the mirror, which is al- 
itiost impossible to obtain using material mirrors in real 
mechanic motion. 

To overcome this difficulty, alternative methods have 
been explored, based on the modulation of the effective 
optical length of the cavity. A scheme based on the modu- 
lation of the skin depth of a semiconductor mirror was pro- 
posed in [B] and is presently being experimentally imple- 
mented [7] . Schemes based on coupling the cavity mode to 
an emitter with time-dependent properties were explored 
in [ilS]. 

The recent observation of a microwave emission from a 
modulated SQUID into a superconducting circuit has been 
considered as a first experimental evidence of DCE [1^ : as 
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the phase of microwave reflection on the SQUID depends 
on the applied magnetic field, a spatially moving mirror 
can be simulated by rapidly varying the magnetic field 
imposed to the SQUID, which has been predicted to lead 
to an appreciable DCE emission [TT| . 

In the present paper, we theoretically explore a strategy 
to observe the DCE in the optical domain by modulating 
in time the refractive index of the medium filling the cav- 
ity. This strategy to modulate the optical length of the 
cavity, originally proposed in [T2HI4) . has recently been 
pushed further in |15j and bears a tight resemblance to 
the concept of time refraction [TB] and to recent efforts 
at using optical analogues for the study of cosmological 
particle creation and black hole evaporation PTIITS] . 

Even though our theoretical developments are very gen- 
eral, here we focus our attention on a specific scheme 
where the optical length of the cavity is modulated in 
time by including in the cavity a slab of x^'^'^ nonlinear 
optical material and pumping it with a periodic train of 
ultrashort optical pulses with an (angular) repetition rate 
r2 = 2-K /T, T being the time separation of neighbouring 
pulses. As a result, the cavity mode (of unperturbed fre- 
quency Wc) experiences an effective refractive index of the 
nonlinear slab which follows the instantaneous pump in- 
tensity ncff = 710 + 5n(t) = no -I- n2lp{t). As usual, the 
DCE emission is maximum when the resonance condition 
fl ~ 2ujr is satisfied. 
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Fig. 1: Upper panel: Proposed experimental scheme for study- 
ing the DCE. The optical cavity resides within a photonic 
nanowire or waveguide and is externally excited by a peri- 
odic laser pulse train that modulates the cavity optical length 
through the nonlinear Kerr effect. DCE photons are excited 
from the vacuum state and transmitted along the waveguide 
which terminate with photon detectors. In addition, the cav- 
ity may be seeded by a optical signal sent down the waveguide 
from one end and detected at the other end. Lower panel: 
simulated pulse train obtained by superimposing the first four 
harmonics of a nanosecond Nd:YaG laser at a fundamental 
wavelength of 1.064 fim. 

To quantitatively describe the DCE emission, a quan- 
tum optical model of the modulated cavity is developed, 
which can be analytically solved and where the details of 
the specific configuration are summarized by a few param- 
eters. Inserting in the model realistic values of the param- 
eters for a pump pulse train resulting from the interference 
of several harmonics of a fundamental laser beam at 1064 
nm, an experimentally sizeable DCE emission in the near- 
IR around a wavelength of 2 /im is predicted. Using the 
same model, we also investigate the response of the cavity 
to a seed beam which is incident on the cavity at a fre- 
quency uJs close to resonance with the cavity mode: in the 
presence of the modulation of the refractive index at f2, the 
stimulated counterpart of the DCE results in the amplifi- 
cation of the seed beam after being transmitted through 
the cavity, and in the appearance of a phase-conjugated 
beam at a frequency 17— due to the scattering of the seed 
beam on the time-dependent refractive index. Our predic- 
tions are finally verified by numerically solving Maxwell's 
equations for a cavity filled with a linear medium with a 
time-dependent refractive index. 

The physical system. — The system we are study- 
ing is based on an optical waveguide cavity as shown in 
the schematic drawing in Fig. [TJ The cavity is embedded 
within an optical waveguide, e.g. as a Bragg grating cav- 



ity. A train of ultrashort, intense laser pulses separated 
by a time interval T is sent externally on to the waveg- 
uide, ideally at 90 deg. The laser pulses, coupled to the 
medium through the x^'^^ optical nonlinearity of the cavity 
material, induce a periodic modulation of the effective re- 
fractive index rioff (^) = no+6n{t) = n()+n2 Ip{t), where hq 
is the refractive index of the unperturbed cavity mode, n2 
is the nonlinear Kerr coefficient (proportional to the x'^' 
nonlinear susceptibility) and Ip{t) is the instantaneous in- 
tensity of the pulse train. 

Using the typical value 712 3 • 10^^® cm^/W for fused 
silica, pulses with a reasonable peak intensity on the or- 
der of 10^'^ W/cm^ will modulate the cavity refractive in- 
dex between hq and uq + 0.003. A pulse train with the 
required intensities can be generated by combining the 
first four harmonics of a nanosecond Nd:YaG laser in a 
similar fashion to pulse train generation techniques al- 
ready demonstrated in the literature [Tallin]- As shown 
in the lower panel of Fig. [Tl this will create a train of 
pulses with sub-femtosecond duration and spatially sep- 
arated by the fundamental wavelength A = 1064 nm 
of the beam, which gives an (angular) repetition rate 
fl = 27rc/A = (27r)0.3 PHz. Less than 10 mJ energy, 
equally distributed between the four harmonics, is suffi- 
cient to guarantee the required intensities. Care should be 
taken that the cavity transverse dimension is kept smaller 
than the distance between each pulse so as to ensure that 
only one pulse at a time modulates the cavity length. Typ- 
ical ridge waveguides, photonic crystal cavities or photonic 
nanowires all satisfy this requirement and are viable solu- 
tions. The time-modulation of the refractive index of the 
cavity will excite real 2 /im photons out of the vacuum 
state by DCE, which can be efficiently collected at the 
output ports of the waveguide and detected with state of 
the art single photon detectors. 

The theoretical modeL — For the sake of simplicity, 
we restrict our description to a single cavity mode with an 
electric field profile 

5,=£,"cos(Me-<"'+^''/"^'. (1) 

fully contained within the fiber. The z-axis is oriented 
along the cavity axis and ctc is the mode waist in the trans- 
verse X and y directions. For a cavity of length Lc filled 
with a material of dielectric constant Ec, the single-photon 
amplitude is equal to £^ ~ ^ ^fnod LcCrl^c- 

The pump beam consists of a train of optical pulses of 
duration t separated by a much longer time interval T 3> 
T, that propagate across the cavity at a speed Vp along 
the direction x, orthogonal to the cavity axis. The peak 
electric field amplitude is E'^. Along y and z, the pump has 
a transverse Gaussian profile of wide waist dp 3> l/fcc, <^c- 
In formulas, this corresponds to an electric field profile 
Ep{r,t) = Ep{r,t) e^'''^p^~"p*^ with carrier frequency ujp, 
carrier wavevector kp, and envelope 

Epir,t) ^ ^ij0g-[-"p(t-"T)lV2^%-(y^-(-.^)/2.^^ (2) 
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Here, the sum over n runs over the pulses forming the 
train, a = VpT is the spatial length of each pulse within 
the fibre material, which is assumed to be non-dispersive. 

The Hamiltonian describing the dynamics of the single 
cavity mode under the nonlinear modulation (Se(r, t) of the 
dielectric constant induced by the pump train of pulses has 
the form |2H] 

(3) 

where Oc and arc the destruction and creation operators 
for photons in the cavity mode. The time dependence of 
the effective coupling constant 



then consists of a train of peaks of height 



Ao = - 
and effective duration 



2 no Lc VpT 



^/al + (T^/vp 



(4) 



(5) 



(6) 



that result from the spatial overlap of the cavity mode 
with the pump beam. Note that the effective duration f 
can be significantly longer than the pulse duration r as 
soon as the waist ctc of the cavity mode exceeds the pulse 
length a = VpT. The peak modulation of the refraction 
index is 

SnP''^'' = na/r"'' = ^ \E"J^- (7) 
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The Fourier transform of A{t) consists of a comb of S- 
peaks spaced by the repetition rate 27r/T and multiplied 
by a broad Gaussian envelope of width proportional to the 
inverse pulse duration f~^. 



i=-oo 



(8) 



In the following we shall neglect the small nonlinear fre- 
quency shift of the cavity mode due to the ajoc terms in ([3]) 
and concentrate our attention on the processes where two 
cavity photons are either created or destroyed by Hamil- 
tonian terms proportional to (aj)^ or d^. These processes 
are strongest when one of the components of A{uj) is close 
to resonance with twice the cavity mode frequency. From 
now on, we shall assume this condition to be approxi- 
mately met for the j-th component at Cl = jfl, and we 
neglect all other components, 

A{t) ^ ^/^e-^'^'/4e-^^*+c.c. = HBe-^''' +c.c. (9) 

This leads to the final form for the isolated cavity Hamil- 
tonian, 



t2 ^ g^Ot ~2 



(10) 



For strong enough modulations 2B > $1/2— Wc, the Hamil- 
tonian ([TU]) for the isolated, lossless cavity predicts that 
the number of cavity photons exponentially grows in time 
at a rate 



r = 2W4B2 



UJc 



(11) 



this exponential amplification of the cavity field will be 
recovered in the section about numerical calculations using 
an ab initio FDTD simulation of the electromagnetic field 
under a time- modulation of the dielectric constant. 

To go beyond this very idealized model and be able to 
describe the steady state of realistic cavities, one has to 
include radiative and non-radiative losses and, possibly, 
the external driving of the cavity by means of incident light 
beams. In standard treatments [26l[27], the Hamiltonian 
for a single cavity mode coupled to external light sources 
is written in the form 



H = Ho+ TlKe [Eincit) + E*^^{t) a,] 



(12) 



where Einc{t) is the amplitude of the coherent laser field 
incident on the cavity and coefficient quantifying 

the coupling of the incident radiation to the cavity mode. 
Losses are then included at the level of the master equation 
for the density operator 

dtp = --^[^1 P] + y {2acpal - pajoc - ajacp}, (13) 

where jc is the total cavity decay rate. 

Analytical predictions. — Standard quantum opti- 
cal techniques such as the input-output formalism or the 
semi-classical Wigner representation of the quantum field 
dc can be used to obtain exact predictions for the most 
significant observables in the steady-state reached by the 
cavity under the interplay of the Casimir modulation, the 
losses, and possibly a coherent incident laser beam pSll?7] . 
In the absence of a coherent drive Ei^c ~ 0, the total pho- 
ton emission rate 



DCE 



j^/4-AB^ + {Lo,-n/2y 



(14) 



vanishes in the absence of any modulation B = and is 
strongest at the DCE resonance, i.e. when the modulation 
frequency f2 is close to twice the cavity frequency ujc- The 
divergence of the emission rate for B — > iJthr with 



1 



UJc 

2 



(15) 



signals the threshold for coherent oscillation in the cavity 
by a mechanism which is the dynamical Casimir analog of 
parametric oscillation. A related parametric oscillation ef- 
fect via DCE was observed in [3T] using a superconducting 
circuit cavity whose electric length is modulated in a fast 
and periodic way by an external magnetic field driving the 
terminating SQUID. 
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Fig. 2: Relative intensity of the transmitted [panel (a)] and 
phase-conjugated [panel (b)] beams as a function of the seed 
frequency ujs- Black solid lines: Q — 2bJc, S/Tc = 0.05. Red 
dashed lines: f2 = 2cjc, -6/70 = 0.15. Blue dotted llines: H — 
2cj, = 27,, Bhc = 0.3. 



In the presence of a coherent seed i?inc(i) = e^*"=* 
at frequency w^, the coherent emission on the opposite 
side of the cavity consists of two monochromatic beams: 
a transmitted beam at the seed frequency cOs and a phase- 
conjugated beam at the frequency — Ws, specular with 
respect to n/2. For a symmetric cavity, the two emerging 
beams have amplitudes 



Et = 



7c/2 



{ijjs - ujc + i^){uJs + - O + i^) + 4^2 



(16) 

K (17) 



Examples of the Wj-dependence of their intensities are 
plotted in Fig. [2] for different values of the modulation 
amplitude B and of the detuning H — 2aJc- 

In the absence of modulation S = 0, the amplitude Et 
of the first component [Fig. ^a)] reduces to standard res- 
onant transmission through the cavity: it is complete on 
resonance uJs = and the resonance peak has a linewidth 
7c. For finite modulations B, the transmitted intensity 
can grow above one, which signals the onset of stimulated 
DCE processes. For increasing values of the modulation 
amplitude B, the resonance peak in the Ws-dependence of 
the transmission becomes sharper and the linewidth of the 
amplification peak tends to zero as the oscillation thresh- 
old is approached 

A similar ajs-dependcncc is apparent in the amplitude 
Epc of the phase-conjugated beam that emerges from the 



-"^ Simple algebraic manipulations show that Et{u}s) and Epc{(^B) 
exhibit a pair of resonant transmission poles as a function of uie . In 
the simplest case Q = 2lUc, the two poles are at ojc — j(7c/2 — 2B) and 
f^c — i(7c/2 -I- 2B). In particular, note how the imaginary part of the 
first pole tends to zero when approaching the parametric threshold 
B 7c/4. 



cavity at the specular frequency — ujg [Fig. [1Kb)]. The 
amplitude of this latter beam vanishes in the absence of 
any modulation, so that its very presence is a consequence 
of a dynamical Casimir mixing process of the incident field 
by the temporal modulation of the refractive index. 

For Cl = 2ll!c, the Ws-dcpendencc is also peaked at cjc and 
shows a linewidth narrowing phenomenon as the threshold 
is approached B — > 7e/4. In the presence of a finite detun- 
ing of the modulation 57 — 2lUc 7^ 0, the transmission peak 
splits at low B into a doublet, one peak being close to Wc, 
the other peak being close to — oJc- As B is increased, 
the peaks move closer to each other and finally merge. As 
the threshold is approached, the splitting transfers into 
the imaginary parts as in the = 2ujc case. 

We conclude this section with some quantitative re- 
marks about the actual intensity of the emitted DCE ra- 
diation. To this purpose, we insert realistic parameters 
into the expression ([Ti|) for the photon fiux emerging from 
the cavity. From the definition of i3 in ([5]), it is immedi- 
ate to see that the result is dominated by the Gaussian 
exp(— fj^f 2/4) and a critical condition to have an appre- 
ciable photon generation is that that f < 2/lOc = Tc/'k: 
the effective duration f of the Gaussian pump pulses in 
the cavity has to be comparable or shorter than the op- 
tical period Tc of the cavity mode. Otherwise, the cavity 
mode sees the modulation of the index as almost adiabatic 
and remains in the ground state with no photons. On this 
basis, it can be advantageous to tune the repetition rate 
close to twice the cavity frequency, ~ 2a;c, so to work 
with the fundamental component of A{uj) for which the 
Gaussian factor is least crucial. 

Moreover, for a given value of 5iiP°^'^ and a fixed r/T, 
the uJc factor in (O suggests it is advantageous to shift the 
cavity frequency uJc (and correspondingly the repetition 
rate Vl) towards shorter wavelengths, e.g. in the visible 
domain. Indeed the number of cavity photons is predicted 
(in the lossless case described by Eq. ((TT|) and considered 
in the numerical calculations of the next section) to grow 
exponentially with B cx lOc- However, there is a trade-off 
between increasing the frequency and the actual possibil- 
ity to achieve the necessary modulation rates. Trains of 
^ 1 fs pulses with a few micron periodicity have been 
demonstrated and sub-fs durations with 1 micron peri- 
odicity should be relatively easily generated as discussed 
above. However, shorter periodicities would require shift- 
ing the spectrum further into the UV region, thus encoun- 
tering issues with material absorption and a general dif- 
ficulty in obtaining high energy pulses e.g. at the fourth 
harmonic of an 800 nm laser pulse. We therefore believe 
that, whilst the laser source scheme adopted here is per- 
fectly reasonable, it docs also seem to be the limit allowed 
by current technology. 

Realistic values for the pump train of pulses can be sum- 
marized as follows: combining the first four harmonics of 
the fundamental wavelength at A = 1 /.tm, a pulse dura- 
tion on the order r = 1 fs can be obtained. Realistic 
values for the cavity parameters can be 5nP°^^ = 0.003, 
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Fig. 3: Numerical simulations for a cavity excited exactly on 
resonance: = = (27r)0.3 PHz, 5n = 0.003. Snapsiiots (a) 
of the electric field inside the cavity for a series of increasing 
times and field intensity versus time (b) for a cavity of length 
Lc = 2 fim. (c) Logarithmic gain versus modulation frequency 
SI for three different cavity lengths, 8 fim (solid line), 10 fj,m 
(dashed line), 16 fim (dotted line). The modulation is always 
taken on resonance with a cavity mode, f7 = 2ujc. 



<7p/ Lc ~ 0.5, 7c = 6 • 10^"'^ s^^ and CTc = 0.5 /iin, which 
leads to the quite optimistic prediction $dce ^ 10^° s^^ 
for the emitted photon flujo. Assuming that the laser 
pulse train lasts for 1 ns, we predict a DCE photon emis- 
sion rate of the order of 10 photons/pulse that holds strong 
promise for detection. And it is worth noting the value of 
B is also not far below the threshold value Bthr for DCE 
parametric oscillation. 

Numerical calculations. — The calculations of the 
previous section were based on a quantum optical descrip- 
tion of the evolution of the field amplitude in a single cav- 
ity mode coupled to a dissipative bath of external radia- 
tion modes. A most significant advantage of this approach 
is the possibility of including non-linearities in the model, 
so to describe the back-action effect of the DCE onto the 
external modulation p2l[23] . 

In this section we explore a different strategy to the- 
oretically study the DCE, based on the FDTD solution 
of the classical Maxwell equations with a time-dependent 
refractive index. Zero-point noise in the initial quantum 
vacuum state is included in a phenomenological way via a 
noisy initial condition. A key advantage of this approach is 
that it is able to follow the system evolution in real time 
and therefore quantify the actual time-scale over which 
the DCE gain can be observed. Another promising point 
is the possibility of extending the calculation to cavities 
of arbitrary geometry, so to fully include effects stemming 
from the multi-mode nature of the field. 



^Note that for a propagation speed iip < c, a pulse duration of 
1 fs corresponds to a length a < 0.33 fim. An important contribution 
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Fig. 4: Numerical simulations of the DCE in a cavity of length 
Lc = 80 fim seeded by noise. The modulation at £1 = {2tt) 0.3 
PHz is resonant with a cavity mode and has Sn — 0.003. (a) 
input spectrum at t=0 ps (red line) and spectra at two different 
times, t = 10 (blue line) and 13 ps (black line), (b) Electric 
field distribution inside the cavity (zoomed in to the first 5 ^m) 
for six different successive times starting at t=0 ps (black line), 
(c) Temporal evolution of the field intensity inside the cavity. 



We start by considering the simplest geometry of a 
spatially homogeneous cavity material of refractive index 
uq — 1.5 enclosed between perfectly reflecting cavity mir- 
rors. These are simulated by imposing perfectly reflecting 
boundary conditions to the field [24]. The cavity index is 
made to sinusoidally oscillate with an amplitude 6n and 
a seed pulse is inserted as an initial condition. Figure [3] 
shows an example for Sn = 0.003 and a Lc = 2 jum long 
cavity with a seed pulse wavelength resonant with the cav- 
ity mode that has a wavelength (in vacuum) Ac = 2 /im 
and excited under the resonance condition D, = 2a;c. In 
Fig. ^a.) we show the increasing electric field amplitude 
for different times and in Fig. |3tb) we show the field in- 
tensity inside the cavity as a function of time. In Fig. EJc) 
we show the logarithmic gain, 
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(18) 



to T then comes from the mode waist 



for three different cavity lengths Lc and varying modula- 
tion frequency f2: as expected from the model, the gain 
does not depend on Lc- On the other hand, as is in- 
creased, a linear increase of the gain, implying an expo- 
nential increase of the field intensity inside the cavity is 
observed, in agreement with the analytical calculations. 

In Fig. |4]we show another simulation in which we insert 
noise as an input condition: this phenomenological way of 
including in a classical FDTD simulation the zero-point 
fluctuations of the fleld in its quantum vacuum state can 
be made rigorous in terms of quantum Langevin equations 
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Fig. 5: Numerical simulations of the DCE in Bragg grating 
cavity, (a) grating index profile and electric field amplitude at 
t=30 ps. The cavity index is modulated at SI = {2n) 0.3 PHz 
and 5n = 0.003. (b) The field intensity evolution with time 
inside the cavity. 



or using the Wigner representations of the quantum field 
P51l27| . Figure m shows the initial spectrum (t=0 ps) and 
the spectrum after 10 and 13 ps (a), and the electric field 
distribution at different times (b - the thin black line shows 
the initial "noise" condition) with Sn = 0.003 and A = 
27rc/r2 = 1 /im. A single peak is clearly amplified at the 
resonance condition Wc = = (27r)0.15 PHz. The fact 
that that a steady exponential growth sets in after just 
^ 10 ps, as seen in Fig. |H[c) is an important indication of 
the accuracy of the analytical single-mode model. 

Finally, in Fig. [S] we show a simulation in which we sub- 
stitute the perfect cavity mirrors used above with Bragg 
grating mirrors: the cavity is thus composed of two 30 
layer Bragg reflectors (sinusoidal refractive index variation 
between 1.45 and 1.55) that enclose a 170 fim long cavity 
that has uq ~ 1.5 and Sri = 0.003. The cavity length was 
fine tuned so as to exhibit a resonance at 2 /im wavelength. 
The cavity and the electric field distribution at the end of 
the simulation are shown in Fig.jSJa). The intensity evolu- 
tion inside the cavity modulated at resonance (A = 1 /im) 
is shown in logarithmic scale in Fig.[5l^b): after ~ 10 ps a 
steady growth regime is reached with a rate of ~ 4 dB/ps. 
The oscillations in the field intensity appear to be related 
to the beating of the modulation with the propagation of 
the pulse back and forth in the cavity with a periodicity 
that depends on the cavity length. With very short cav- 
ities this effect is not observable, but longer cavities lead 
to this characteristic beating with a periodic cavity mod- 
ulation. Similar behavior is obtained with a wide range of 
input parameters (e.g. cavity lengths, Sn). 

Conclusions. — In conclusion, we have proposed and 
quantitatively characterized an experimental scheme for 
measuring the dynamical Casimir effect: the effective op- 
tical length of a cavity is varied in time on a PHz time 
scale by means of a train of ultrashort pulses that mod- 
ulate the effective refractive index of the cavity material. 
The predicted rate of photon emission for realistic pulse 
train and cavity parameters appears promising in view of 
experimental studies. The over-arching idea of laser pulse- 
induced cavity modulation is very general and may be ex- 
tended to other settings such as acoustic cavities modu- 
lated by laser-pulse trains 
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